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HEAT TRANS PER OP AIRFOILS AND PLATES* 
By Otto Sei"bert 

The few available test data on the heat dissipation 
of wholly or partly heated airfoil models are compared 
with the corresponding data for the flat plate as obtained 
"by an extension of Prandtl's momentum theory, with differ- 
entiation "between laminar and turbulent boundary layer 
and transitional region between both, the extent and ap- 
pearance of which depend iipon certain critical factors. 
The satisfactory agreement obtained justifies far— reach- 
ing conclusions in respect to other profile forms and 
arrangements of heated surface areas. The temperature 
relationship of the material quantities in its effect on 
the heat dissipation is discussed as far as is possible 
at the present state of research, and it is shown that 
the profile drag of heated wing surfaces can increase or 
decrease with the temperature increase depending upon the 
momentarily existent structure of the boundary layer. 



INTRODUCTION 



Economical high— speed flight requires aircraft of 
high aerodynamic quality - that is, aircraft with air 
f r ame des igned for minimum flight drag and all parasite 
areas, especially the radiators- reduced to a minimum. 
In this endeavor attempts have been made in this country 
as well as abroad to xxtilize exposed outside surfaces, 
preferably those of the wings , for heat transfer. Such 
radiators are hereinafter called "skin radiators." 

The heat output of the exposed outside surface of a 
skin radiator and its distribution across the surface are 
hardly amenable to prediction from the differential equa- 
tions for flow and heat exchange, as these equations are 
not integrable in such generalization. However, the 
moment theory and its refinement according to Prandtl af— 
fords at lea st very satisfactory practical average values 

* M ¥armeubertragung von Prtfilen und Flatten. " Jahrbuch 
1938 der deutschen Luf tf ahr tf ors chung , pp. II 345 - II 256. 
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of the heat transfer coefficients. The only stipulation 
for its validity is that skin friction and heat transfer 
"begin and stop simultaneously. In the practical case 
that some portion of the outer surface of an airplane, 
whether on the fuselage, engine nacelles, or wings are 
designed as skin radiator, the heat effect usualls* - begins 
after the friction; hence a temperature profile forms 
progressively in the already existent and still growing 
frictional boundary layer which unfortunately does not 
lend itself to representation in closed form, although it 
is mathematically defined by the equationsof motion and 
heat exchange. The correlation between the Nusselt num- 
ber Hu and the Reynolds number Re for all possible 
flight conditions must therefore be secured by measure- 
ments on models. The subsequent study deals primarily 
with the processes on wing radiators — that is, radiators 
mounted in the wings. In this connection the data on 
the heat dissipation of flat plates are invaluable; in- 
deed, they may even save special measurements on airfoils 
(see sec. V) if no unusual claims on accuracy are involved. 



II. AVAILABLE TEST DATA 



The most important publications on investigations 
of this nature, all of which deal with wing radiators, are: 

a) R..&M. Report Ho . 1311 ( 1927) 

Vind— Tunnel Tests on G-loster and Supermarine 
Wing Radiators (Supermarine airfoil section 
RAF 30, chord 15 ft = 1524 mm, model scale 
1:1, span = tunnel diameter, of which 34 in.= 
8 76 mm were heated, total upper and lower 
surface heatable with water separately or 
collect ively ) . 

b) R.&Ii. Report So. 1163 ( 1928 ) 

On the Convection of Heat from the Surface of 
an Airfoil in a Wind Current (model of RAP 
section 26, 5 in. = 152 mm chord between end 
plates at 3 in. - 203 mm spacing, with 29 span- 
wise platinum strips over the total surface). 

c) R.&Ii. Report Ho. 1326 ( 1930) 

Wind— Tunnel Experiments on Steam Condensing 

Radiators (nose st earn— heat ed , RAE 30 section, 
2 ft = 610 mm chord, 7 ft = 2134 mm span, 
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the radiator occupies 5 ft = 1524 mm; 25.1 
percent of the arc length of the upper sur- 
face and 15.2 percent of the lower surface 
or 24 percent and 14 percent, respectively, 
of the chord (of the symmetrical) profile 
heata"ble with saturated vapor at low posi'— 
t ive pr es sur e ) . 

d) E.AM. Report No. 1481 ( 1932) 

Estimation of Wing Surface Area for Evapora- 
tive Cooling (comparison of R.&M. Reports 
Hob. 1311 and 1326). 

In figure 1 the Nusselt number Nu = f (He) is shown 
according to R.&M. Reports Nos. 1163, 1311, 1326, and 
1481 for different angles of attack B. 



wher e 

c 



Nu = ~— - Nusselt number 

o 



Re - ■—— = Reynolds number 

D airfoil circumference 
v flying speed 

v kinematic viscosity referred to state of air in un- 
disturbed zone 

\ 0 heat conductivity of air at surface portion. 



t 0 ~ 40° on RAP 26, 
t 0 ~ 9 3° on RAF 3D; 

Air temperature assumed at 20°. 



III. CONCLUSIONS 



a) The heat emissivity of the surface portions in 

proximity of the leading edge is greater than 
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at any other part of the airfoil. It naturally 
increases with Re, "but varies so much lees 
with the angle of attack 3 as the area ex- 
tends less on the upper and lower surface. 

To) On the upper and lower surface the heat dissi- 
pation decreases along the wing trailing edge 
for given Be; "but in the transitional zone 
from laminar to turbulent "boundary layer it 
can increase once more and then decrease, as 
exemplified in the test data of R.& M. Ho . 1163. 

c) By ascending He and g the heat dissipation 

on the upper surface increases, the function 
Hu .= f(Ee) resembles on the whole a curve at 
at least as steep as that for the nose. 

d) The heat dissipation on the lower surface of the 

wing scarcely varies over the practical range 
of positive angles of attack; it increases at 
angles of attack corresponding to zero or 
negative lift. Depending upon the absolute 
magnitude of Re (and other factors, see 
pt. IT) the curve can be as flat as that of 
the nose or about as steep as or steeper than 
that of the upper surface. 

While this dependence of the heat dissipation is 
in general quite comprehensible, the different slopes of 
the curves for the heat dissipation of the surface por- 
tions cannot be summarily explained. These and other re- 
lations are discussed in part IV in the light of the en- 
larged momentum theory as applied to the flat plate. The 
comparison of these theoretical relations with the results 
obtained on the RAP 26 and RAF 30 sections (fig. l) be- 
comes part i cular ly simple. The numerical agreement is 
also very satisfactory if allowance is made for the dif- 
ference in the flow velocity on the upper and lower wing 
surface areas. The arguments can likewise be applied 
quite satisfactorily to heated leading edges by a minor 
conversion. Reliable heat dissipation data on flat plates 
are available (reference l) but they do not lend them- 
selves to strict comparison with the ideal case, because 
the test arrangement causes the heat dissipation to start 
after the skin friction. In any event they enable to 
support the theory and, in addition, afford some insight 
into the expected displacements when only parts of the 
wing surface area are heated. Very insuff icient, however, 



NACA Technical Memorandum Ho. 1044 



5 



are the data available on the effect cf the temperature 
increase of the heated portion or the- total surface, re- 
spectively, or better, of the dimensionless surface t em— 

perature — (see pt . IV) on the heat dissipation and, 

above all, on the drag of the total area. This might be 
regarded as pure skin friction on the flat plate, while 
on the airfoil the heating itself can influence the form 
drag (pt. VII). 

IV. HEAT DISSIPATION OTP FLAT PLATES IH THEORY AND TEST 



On real flat plates— that is, those of finite thick- 
ness- the boundary layer is theoretically formed some- 
what differently than on ideal plates — that is, on infinitely 
thin plates in symmetrical flow parallel to its plane. 
Postulating complete freedom from friction the ideal plate 
does not affect the flow at all, while the real plate out- 
wardly displaces the individual streamlines, depending 
upon the shape of the leading edge, and so causes increases 
of speed relative to the undisturbed flow. 

Taking into account the skin friction, the boundary 
layer on the ideal plate is produced by the decelerating 
action of the surface particles on the flowing fluid which 
then continues throtigh it. Hence in two-dimensional flow 
the velocity profile of a section at right angles to the 
flow direction at any distance from the plate leading 
edge must (as on a pipe) bo unequivocally defined by the 
Reynolds number (reference 2) computed with this depth x 
and the undi stur bed velocity v and the kinematic vis- 
cosity V i 

Rex = ^f- (1) 

regardless of whether and how far the plate extends down- 
stream. On the real plate this skin friction is superposed 
by the displacement process associated with the shape of 
the body (especially its thickness, over— all length, round— 
ings ) which very likely can also become effective upstream. 

The combined action of this pr of i 1 e effect and the friction 

favors the creation of the so— called form resistances which 
cannot occur on the ideal plate. 
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The velocity field on the real plate must accordingly 
have a somewhat different aspect, which is, that the dis- 
crepancies from the theory are especially great at the 
"beginning of the plate. The effect of the leading edge on 
the heat transfer is obviously similar. The foregoing 
logically applies also to airfails with and without angle 
of setting, where an even somewhat poorer agreement is to 
be expected. 

On ideal plates and on real thin plates* (with sharp- 
ened edge and int er f er en ce— f r ee inflow) a laminar flow 
forms first in the "boundary layer. Turbulence does not 
come into being until after traversing a certain entrance 
length x cr which can be computed from the so— called 

critical Reynolds number. 

Rex cr =^l (2) 

Por x > x cr the laminar boundary layer becomes fully 
turbulent in a transitional region; Rex cr according to 
Ten Bosch (reference 2, pp. 139—140), has been measured up 

to 5X10 for very small initial disturbances, but may drop 

to 10 5 in the presence of major disturbances, the transi- 
tional region can extend from Rex = 10 5 to 2X10 7 for 
sharpened edge and steadied flow. With the one— sided sur- 
face C of the plate, the drag coefficient Cf and the 

air density p the drag of the ideal plate V for each 
side is 

V = OfO^E (3) 

The amount of heat dissipation from one side of the plate 
at increase of temperature * 

Q = a mK 0* (4) 



is, according to simple momentum theory, and as employed 
by Yon Karman and Latzko (references 3 and 4), to the 

*Thus Ten Bosch (reference 2) indicates on page 148 that 
heat exchange surface of streamline form should have less 
than 33X100/600 = 5.5 percent profile thickness, since 
only from this slenderness ratio on the friction forms the 
principal portion of the total drag. (See also pp.139— 140.) 
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plate drag as the heat content of some gas mass in the 
nucleus to its quantify of motion: 

A = LEI ) ( g c pm * ) _ g c pm» / - \ 

¥ - Tpvjv v k 

Equations (3), (4), and (5) afford 

A <V 

<%K = i~ = i0 f c pm vf (6) 



or the mean Uusselt number with the plate chord x in- 
troduced in equation (l): 

= I- C f (*5jHa)(i££)-0 f i- Pr He, ( 7 



A similar line of reasoning was originally followed by 0. 
Reynolds (reference 5) to the flow in circular pipes. But 
here, as on the flat plate, this simple theorem enables a 
satisfactory reproduction of the actual amount of heat 
transfer only in the specific case where the Prandtl number 

gn c p 
Pr = - x - 

is equal to unity. In all practical cases for which 
Pr ^ 1; hence for air also* divergencies are found which 
Prandtl attempts to comply with by assuming a special 
gas fi lm between surface and actual boundary layer within 
which the gas velocity increases in proportion to the wall 
distance. Its thickness and velocity distribution existing 
in it are defined by the proportionality factors and cp , 
to be computed from experiments (reference 2 ) / Hence, ac- 
cording to the enlarged Prandtl equation: 

_ \ Cf c pm v V 
a m = * = cp(l-PrJ (8) 



*Pr~0.725 for -air of surrounding temperature; Pr varies 
only slightly with temperature. 
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ITux = 



1_ 

2 



Cf 



\|/ - cp (1 



c td m ^vx 



Pr g ) 



= 0f- 



•^Pr 
2 r 



\|/- cp(l - 



Pr g ) 



Hex 



(9) 



with P*"g indicating the Prandtl number for the " 'bound- 
ar y laye r t en roer at ur e . " Prandtl gives the values \|/ and 
cp only for the round pipe. In the case of 



heating of the fluid = 1 cp = 1.4 



Pr -o.ia 5 Re -o.i 



cooling of the fluid = 1 



and in "both cases 



Re = 



9 

vd 
v 



=1.12 Pr 



-0.185 



Be 



-o.i 



d = pipe diameter 



— 0.185 



0.075 



Ten Bosch comouted \|/ = 0.89 and cp = 1.4 Pr Hex 
from experiments for the flat plate. The heat transfer 
factor is further dependent upon the difference 9 be- 
tween surface temperature and gas temperature. Ten Bosch 
assumes a turbule n ce v ar iat ion ; \ as factor in the drag 
coefficients. Lacking more precise data he gives quantity 

1 as function of the increase of temperature 9. Figure 

2 holds true for flow in round pipes. Ten Bosch further 
suspects a probably minor dependence on Pr . For a more 
exact exploration of these important interconnections the 
author recommends in place of 9 two variables: the abso- 
lute temperature T of the gas and the ratio ^g/^ — that 

is, the " dim ens i onl es s 11 boundary layer temperature - where 
Tg is the mean absolute temperature of the boundary layer. 

If, as in the present case, only one gas is involved, the 
absolute wall temperature T w may be substituted for 

Pr discounted; 9 



as determining quantity, and 
tive for heating of the air. 



is po 



g 

, i- 



Since the dependence l = f(9) for the flat plate is 
not known, figure 2 will serve a« a basis within the frame- 
work of the present report. 



In relation to equation (7), equation (9) 
the imxiro vem ent factor 



contains 



V = 



ii - 9 (1 - Pr g ) 



(10) 



which for constant Pr and Pr g can be represented as 

(slightly variable) power function of Hex ,(cf. equation 
(18)). 
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The equations for Of', oc^q » and Nux - in the three 
zones: 

a) In Laminar Boundary Layer (Subscript L): 
For isothermic flow 

c fLis = 1-328 Hex -0 " 5 
for any increase of plate temperature 

C fL = 1.328 i L Rex~°' 5 ( 1 1 ) 

with equation (ll) and the simplification through Pr g = 1 
in equation (8), we get 

1.328 % ^ Rex -0 ' 5 vc pm V 0.664 . , p r , 

«BL^r- f = ; * L X "T" ^ Rex 

against the approximation according to theoretical solu- 
tion conformably to E. Pohlhausen (reference 2, p. 144) 

0 x B m 

Tor the completely free plate, =1; t is unknown in 

Bm 3 

the first equation, and, although the term 7 Pr in the 
second equation represents an approximation the 'relia- 
bility of which is limited to the range of Pr = 1000, it 
is nevertheless used for air. Besides, both equations 
become identical for \|/ = 1 , |_ = 1 and Pr = 1. Ex— 

Jul 

periment alone can indicate which relation comes closest 
to reality. A variation of the "turbulence" in the laminar 
boundary layer owing to finite increase of temperature is 
to be disregarded. It might be assumed that I becomes 
less than 1 for heating (positive 8), as in turbulent flow 
in pipes, and greater than 1 for cooling. 

Henceforth: 

Eux L = -rj- C f -/Pr Rex 0,5 , Pr - 0.725, i- VVr = 0.449 ( 12) 
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"b) Transitional Region inJBoundary layer 
( Subs cr ipt U ) 

According to an approximation "by Prandtl the drag 
coefficient of the t o/^t^l plate for Rex "between 5X10 
and 2 X 10 7 — % allows for the increase of temperature — is 

The figure 1700 is merely the mean value for a range with- 
in which the true value can fluctuate, depending upon the 
initial disturbances.* Although all important conclusions 
could be drawn from the variable slope of this curve from 
its lower limit of validity (Rex -5X10 ), an approxima- 
tion by Ten Bosch is introduced which, 'while not quite as 
close, has the advantage of bringing out the effect of 
Rex cl . on the drag coefficient and naturally also on the 
heat dissipation. 

The total frictional drag is visualized as being com*— 
posed of the forward plate portion L on which the bound- 
ary layex* is laminar, and of the rest of the plate portion 
U , along which the boundary layer transforms, whereby the 
latter portion might be put at 

0 f u = 0.003 ( 14) 

independent of Rex. With equation (10) we get in this 
instance: 

|c f §Pr 

Nux u = ^-ryxr^FrTT Sex = c ™ Y r Pr Re * (15) 



c) Turbulent Boundary Layer (Subscript T): 
lor Rex > 2 X 10 7 



*T'en Bosch used an older form giving almost the same values 
as equation (13), the first term of which reads 0.072Rex _O " 
lor further information on the drag coefficients of flat 
plates at large Reynolds numbers see reference 6. 
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O fg i = 0.072 \ Sex °' S ( 16) 

By observance of the above values for ^ and cp and 
approximation Pr~Pr g , the improvement factor according 

to equation (10) 'becomes 



V-1.36 for Rex = 5X10 5 

V -1.30 for Rex = 10 7 

or, represented as power function: 

Y = = 1.638 Rex" 0 ' 0143 (18) 

vR ex 

This relation needs to be checked by accurate tests and 
improved accordingly. 

Jigure 3 shows the drag coefficients C f according 
to equations (ll), (13), (14), and (16) and the mean 
ITusselt number Nux over the range Rex = 2,5 X10 4 to 
1.5 X 10 for the unheated plate. Posting the respective 
term for C f and of V in Hux according to equation 
(18) (again for 9=0) affords 
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a) for the laminar 'boundary layer according to equation (12) 

SFux L = 1.328 X 0.449 Sex 0 * 5 = 0.59? Hex 0 ' 5 
"b ) for the transitional region according to equation ( 15 ) : 

HTuxff = 0.003 X 1.638X0.3625 Hex 0,9857 

= 0.00178 Hex 0,9857 

or, according to Prandtl, 

w - 0.455X1. 638X0.3625 Rex 0,985 7 

Xpr (log Hex) 375 " 8 

- 1700 X 1.638 X 0.3625 Hex' 0 - 0143 

O 9 8 5 7 

Fuxp- = 0.27 ----- 1010 Rex" 0 - 0143 

rr ( log Hex) 0-858 

c) for the turbulent boundary layer according to equa- 
tion (17) 

ffuxj = 0.072 X 1.638 X 0.3625 Rex 0-7857 

= 0.0428 Rex 0 - 7 857 

according to equation (?): 

iTux K = 0. 072 X 0. 3525 Rex 0,8 = 0.0261 Rex 0,8 



I or easier plotting of the individual theoretical 
j ines figure 3 contains several numerical values follow- 
ing from these equations. 

The reason for the choice of the specific case of in- 
finitely snail increase of temperature on the plate surface 
Y as *he aT °idance of uncertainty in the choice of % which 
in this instance "become equal to 1. The drag coefficient 
of plates with finite increaae of temperature is theoreti- 
cally analyzed in part 711. For Rex <Hex cr 

Cf ~ Ofi; - Nuxt 

For Rox>Rcx cr up to the end of the transitional region 
(at Rex -2X10 7 ), the moan values arc: 

Cf l[Rcx cr l Box cr + C f -j (Rex - Rox cr ) 

c f - " 

Rex 

= (Cfl[Rox cr l- Qftf) + Gfjj 

^u.x L fHexn-rl Kox cr + Nuxfj (Rcx - Rex cr ) 
J>i UX = — 



Rex 

= (Hux L ^ HoXcr1 - Fuxy) + Nuxff 



R cx cr 



UACA Technical Memorandum No. 1044 



13 




ratio 



is at the same time re- 



The two "cases 



Rex-,. = 10 5 and fiex_ = 5X10 5 



are dif f er-ent iat ed . 

On plartes of such total depth x that Rex ranges 
"between 10 5 and ~2X10 6 , Cf and Hux values situated 

in the shaded region should he possible to "be measured 
(up to. ±40 percent referred to some mean value). (See-ref- 
erence 7.) For still greater plate depths the uncer— 
tainty decreases then rapidly (figure 3). Prandtl's 
formula gives a curve similar to that of the linear com- 
position under the assumption of especially great freedom 
from disturhance. 

The heat dissipation (and the temperature and velocity 
fields) of a flat plate have teen treated in detail "by 
Elias (reference l). The plate was. 500 millimeters deep, 
250 millimeters wide and 29 millimeters thick with a 
wooden frame around it. The strip facing the flow was 
x x - 100 millimeters deep and carefully tapered for the 
purpose of minimizing the initial d i s tur "ban ce s . Referring 
the heat volume Q, , dissipated from "both sides up to a 
certain depth x, to the unit span, he plotted the (named) 

term against Rex (or Pex). This is then compared 

with the theoretical value — CfCp m Y following from 

equation (6). The displacement of the start of heat dissi- 
pation-relative to that of friction he attempted to ac- 
count for "by measuring the depth of the respective section, 
first, from the beginning of the thermal reaction as x^ 
and then from the "beginning of the hydrodynamic reaction 
as x^. But neither of. the representations is really able 
to give a satisfactory reproduction of the process. 

On the other hand, the reasoning on the basis of the 
following is very informative despite the marked scatter: 
if the plate had been heatable up to the leading edge it 
would naturally have dissipate.d more heat by reason of its 
larger area and the added existence of a region of strong 
transformation of the velocity profile and hence of best 
heat dissipation. In the dimens ionless representation 
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the upper limit value Nux max theoretically corresponds 

to the curves Nux according to equation ( 19 ) . The ac- 
tually heated portion would, of course, have then "been 
a"ble to dissipate less heat than in the experiment^ namely, 



Eux 



Nux 



•mm 



(x=x) * - ^(X^j) 



lux 



( x-x) 



= Kus (x=x) - Nux (x= Xl ) 
Hex(. x=x ) 

-llUX/ \ J— 

K * Rex (x=x) ' 



> 



(20) 



fhe measured heat or its corresponding llux m . must lie be- 
tween t v esp two limits. 

Jior comparison with the test data of Uux m , the value 
Uux T is used in figure 4. 



For a constant x x /x equation (20) becomes 

0.7857 / x NO. 7857 

Nux min = 0.0428 Hex/ _ , - 0.0428 ( ^£ ei K 4 

k x " x ' \ x /(x=x) 



ill 



0.0428 Rex 



0.7857 



k (x = x) L 
= Hux T(x = x) [1 - (%") 



x M . 7857 



> (21) 



1 .7 857 



The factor w -that is, the ratio of the actual mean 
Nusselt numher referred to- that of the total depth 



w = 



Hux m 



(22) 



1 .7 85 7 



ITuxj 

therefore ranges between 1 — 

and 1. .The 

parameters are x^x = l/ 2 , l/3, 1/4, and l/5, each of 

which is applicable to one plate each with a "cold" forward 
and a hot" residual portion amounting to 2, 3 4 and 5 
times that of the front portion each in a different air- 
speed flow v (between 10 and 35 m/s). It is readilv 
seen that, * 

(a) WdfcktascttejrbiLotL-Qf oth^oraneauof- -small Reynolds 
numbers the lTu% are located above the 
Nux K , as anticipated. 
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(t>) The experimental curves (and the w factor^) 

are located between the theoretical limiting 
• values . 

( c) They are, admittedly, not quite so systematically 

staggered according to parameter as equation 
(21) postulates. 

(d) With increasing Hex they approach, contrary to 

equation (21), the curve of Nuxqi - that is, 

the factor w depends on Rex owing to the 
dissimilarity of correlated velocity and 
temperature fields. 

(e) The marked steepness of the actual Uux curves 

from the beginning indicates that the critical 
Reynolds number according to equation (2) is 
close to the lower limit 10 . 

In this comparison, especially for point ( d), it 
shoxild be noted that the function Uuxrp is not as yet 

completely defined as to magnitude and steepness. System- 
atic discrepancies in figure 4, which would have to be in- 
terpreted as temperature effect, are, of course, not dis- 
cernible by reason of the smallness of the chosen 0. 



Y. APPLICATION TO AIRFOILS 



These studies on the heat dissipation of flat plates 
can be applied to airfoils. It is advisable to treat the 
upper and lower surf aces s eparat ely by introducing as char- 
acteristic length x the measured distance around the 
airfoil from the stagnation point, the "start" of the 
boundary layer, up t o the end of the heated area. 

By reason of the variation of the airspeed along the 
aiffoil circumference as a result of the angle of attack 
and of the finite airfoil thickness, individual mean ve- 
locities v x for upper and lower surfaces must be intro- 
duced, which follow from the undisturbed velocity, say, 
with the aid of the mean differences Ap of the static 
pressures on both sides of the airfoil in ratio to the 

dynamic pressure q = — — £. The mean Ap values are 
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obtainable from pressure distribution curves. Assuming 
that the acceleration or deceleration follows Bernoulli's 
equation, we get (for constant air density p): 

2 



— is , of course, posted with its algebraic value (nega— 

q. 

tive for the lower surface at positive lift). 

For simplification the slight travel of the stagna- 
tion point with the angle of attack and the identical 
absolute variations of the lengths x are discounted. 
The pressure distribution need not be known, if satisfied 
with the admittedly very rough approximation that the 
velocity on the wing upper surface is greater by the 
amount Av meters per second, and on the lower surface 
smaller by the same amount than the undisturbed flow 
velocity v. Actually, local velocities amounting to 
multiples of the undisturbed flow prevail on the forward 
portion of the upper surface especially on thick and high 
incidence airfoils, while the decelerations on the lower 
surface are substantially less. On the other hand, 



v^p 



= Ap 



Iz i + Ap 

v ' q 



> 



(23) 



J 



(v + 



Av) 2 



- v2 11 



+ I V3 



= SAp 



4 v A v 



- 2a 



where EAp denotes the mean pressure difference between 



upper and lower surfaces, 
of attack |3 we get 



A p 

A v 

v 



With v i n g area 
= l5lAp cos p 



3P and angle 



-a 



cos p 



v • 



_P_ 

2 



-a 



4 cos 
1 ± 



4 cos p 



(24) 
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Tor normal angles of attack cos 
v x /v are shown plotted against 



0 = 1. The functions 
c a in figure 5. 



The Reynolds number used in the following is obtained 
from the Re employed in figure 1, by correction of the 
velocity from • v to v x according to equation (23) or 

approximation (24) and replacement of profile circumfer- 
ence D by chord x measured on the circumference: 



Rex = Re — — 



x 
D 



(25) 



The conversion of the Nusselt number is simpler 



lux = Nu — 
D 



(26) 



In figure 6 the association of figures 1 and 3 with the 
factors of equations (25) and (26) for the flat plate and 
the two airfoils is compared. It is to be noted that the 
transition from v to v x shifts the points horizontally, 

the transition from D t o x in the direction 1:1. 



On RAF section 26 and RAF section 30 x/D is assumed 
enual to l/2 for upper and lower "urface^. The values for 
the nose on both airfoils are plotted on the assumption 
that the undisturbpd velocity v as located b^twepn the 
two v x is more suitable; x/D is avr raged at 0.20 for 
both sides of thf wing. The curvps in figure 6 cover the 
tpst range of figure 1 from beginning to end. 

In the comparison it is to be noted that a different 
t °mpf rature t 0 or increase of temperature 6 is applied 
to ~ach of the three case=: 





to 


e 


RAF 26 . . 


-40° 


-20° 


RAF 30 . . 


-93° 


-73° 


Flat plate 


Air tempprature 


0° 



The fully heated models manifest the following: 
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a) All straight lines lie practically within the 

zone bounded by the theoretical curves Nuxjj , 
Huxy , Euxrp applicable to the flat plate. 

b) On both airfoils the straights representing the 

heat dissipation of the upper surface at dif- 
ferent angles of attack are packed closer 
together, those for the lower surface a little 
farther apar t . 

c) In conformity with the "narrowness" of the zone 

for large Sex on the (symmetrical) B.AJ 30 sec- 
tion the systems of curves for both sides of 
the wing fall especially close together at 
equal angle of attack. 

d) The slope of the separate straights for the EAF 

26 section (with exception of those for UP— 0.9°) 
rises fairly uniformly with increasing distance 
from Nuxj, and ultimately approaches that of 

ITuxj . 

c) The laminar boundary layer persists for some time 
on the lower surface, but the curve for the 
Rex factors above the upper limit of the test 
range for RAP 26 section appears to become 
steeper and, as theoretically stipulated, hugs 
the straight ITux^f until it deflects toward 

ITuxj in the range of the factors of the RA1P 30. 

Heated wing leading edge: 

f) Here the accord with theory is not so good as already 
indicated. The deciding velocity is nevertheless 
very much dependent on the angle of attack as 
reflected by the fairly great distance between 
the straights on both airfoils corresponding to 
the lift coefficients ~0.24 and 0.64. Added 
to this is the effect of especially high in- 
creases cf velocity in this region. The heat 
dissipation of the nose portion on the upper and 
lower surfacesof the R.A3T 26 section could also 
be computed separately from the quoted data (pt. 
II) and evaluated in the discussed manner, where- 
by the pressure distribution might afford even 
I'better" velocities v x and perhaps even better 

arc lengths x by reference to the respective 
stagnation point. 
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Jron d) it can "be concluded that the boundary layer 
on the wing upper surface changes sooner and "becomes fully 
turbulent quicker than expected according to the value of 
Rex computed with v + Av . Undoubtedly this is accom- 
plished by the cited increases of velocity to which larger 
Reynolds numbers also correspond. The total surface area 
of the RAP 26 model is ''rough" by reason of the glued— on 
platinum strips and it is surprising that the boundary 
layer remains laminar over the greater part of the lower 
surface, as indicated quite plainly by the location of the 
respective straight lines. However, the more scattered 
test values for Ug — 0.9° allow it to be inferred that 
the laminar boundary layer is no longer st'able even at low 
angles of attack and lift coefficients, respectively — that 
is, even for very snail differences relative to the flow 
velocity. According to the foregoing it may be expected 
that in the range of Rex << 10 the laminar boundary 
layer on the upper surface also extends over c ont inuotis ly 
greater depths with decreasing Rex, so that the curves 
Sua become progressively flatter (reference 2, p. 153). 
The S— shape course is naturally less pronounced on the 
curves for the wing upper surface than on those for the 
lower, because they are of themselves quite close to the 
asymptote Nuxrp for maximum Rex. 



71. CHOICE OF MATERIAL QUANTITIES 

The sttidies described under part Ila) to d) were made 
on rectangular model 'airfoils (tip to scale of great 

aspect ratio, only a small portion of the span being used 
as heating surface, or on such between end plates (two- 
dimensional flow closely approximated in both cases). The 
airfoil circumference D serving as characteristic length, 
was little more than twice the wing chord. 

The material values A and ri involved in the din en- 
sionless factor_s_ are related to the temperature (about pro- 
portional to «/t). In none of the practically important 
cases of heat— transfer (flow in pipe, flow around pipe, flat 
plate) has there, up to now, been successfully determined a 
"decisive mean temperature" dependent on surface tempera- 
ture, gas temperature, Reynolds number (and Prandtl number) 
to which the material values were to be referred. In R.&M. 
Reports Uos. 1326 and 1481 the Reynolds number Re is com- 
puted with the kinematic viscosity v for outside air 
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temperature t. In Itfu, X is introduced with the value 
corresponding to that of the surface temperature t 0) 

since the processes in the "boundary layer are justly held 
as the essential. -When, furthermore, Ten Bosch recommends 
the kinematic viscosity v = r\ / p corresponding to the 
surface temperature for the plate in forced flow even in 
Rex, since as proved "by experiment, it was "best "suitable" 
for the plate in free flow, this procedure should answer 
the purpose for airfoils as well. It would involve the 
conversion of the "aerodynamic" Rex to the "thermic" Hex 0 




Rex 0 = Sex --v~ ) = Hex 



(27) 



The conversion factor U 



Hex 



H exr 



is simply a 



T o^o 

function of the two temperatures t and t 0 (if this is 

variable a mean valtie substitutes). To avoid confusion 
between Hex and Hex 0 it would be preferable to use the 
respective Pe'clet number Pex 0 = Pr 0 Rex 0 in place of Rex 



the Prandtl ntimber Pr nay be averaged at ~0.725. 



o, 



The 



conversion factor Ur e X — > P e x 0 =0.725 



Tt) 



also depends 



T 0 Tlo 

on t and to only. Accordingly, if the temperature ef- 
fect is taken into account, all straight lines of figure 6 
relating to the airfoils shift toward the left, while Rex 
remains = Rex 0 for the flat plate. But the Fuss el t number 
for the flat plate must be multiplied by the 



factor t 



still unknown 
1 for the case 
flow through 



Assuming that it is smaller than 
of a heat dissipating surface similar to the 
a pipe, the curves must shift downward with rising plate 
temperature and the close agreement of figure 6 is complete- 
ly lost . 



It would afford for the 



RAP 26 



to 
t 



40' 
2C ( 



U 



Rex 



and for the 



u Ro: 



_ 273 +20 1.86 

pr, T = 0.725 X = 0.647 

* cx 273 + 40 1.95 

RAP 30 t 0 = 93° 
t = 20° 
■Pex 0 = 0.725 2ZL±_20 x l i 86 
273 + 93 2.18 



= 0.49 5 
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Accepting, in the absence of other data, the correlation 
of i with the increase of temperature of the pipe for 
the flat plate, the surface temperature becomes 

t 0 = 40° 0 = 20° | = 0.92; t 0 = 93° 9 = 40° \ = 0.88 

Uficx-^Pexo = Pr = 0.725 

The result is illustrated in figure 7. The percentage of 
discrepancies "between the experimental and theoretical data 
is readily apparent. They amount to more than 60 percent 
Con the RAP 30 section). Accordingly, it is very question- 
able whether Ten Bosch's suggested method actually repro- 
duces the true conditions; in any event it would not ex- 
plain what effects the theoretical values for the heat 
dissipation on the recorded, are able to increase. It fur- 
thermore conflicts with the view that the formation of 
the major portion of the boundary layer proceeds the same 
as in isothormic flow, the absolute temperature of which 
is, however, T 0 rather than T. The effect is therefore 
much too great. 

In any case the arguments seem to indicate that the 
temperature effect for every one of the three boundary 
layer structures involved is of a different nature. Ex- 
periments in this zone will require care to assure laminar 
boundary layer over the entire depth of the airfoil or 
plate, and then to secure a turbulent boundary layer on 
the major part of the surface by using surfaces of maxi- 
mum depth (with purposely great initial disturbance (ref- 
erence 7)). In both cases the material values of the un- 
disturbed flow will be introduced in the factors for the 
evaluation (for X. alsot). The plotting of the curves 
Uux = f(3.ex) with the, "temperature effect" as parameter 
affords a formula represented in the form (T 0 /T) n if 

limited to a gas (air), whereby the exponent n-^ for 

laminar boundary layer need not be equal to the exponent 
nj for fully turbulent boundary layer. Both will be very 

much smaller than 1. It might even be possible to split 
these powers wholly or in part in factors and to correlate 
them with the material values in such a way that these 
could be coordinated with the af or e— ment ioned controlling ' 
mean temperature, .as has been accomplished for horizon- 
tal wires and pipes to a certain extent (reference 8). 
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Form resistance does not contribute directly to heat 
diffusion, although very high heat transfer has "been as- 
certained on the l)ac.: of pipes in oblique flow especially 
at large Reynolds numbers, which are traceable to such 
vortices. Up to the present no theoretical formula has 
been published for this type of heat transfer (reference 2, 
p. 150). The induced drag on the wing is, in this sense, 
equally to be treated as form drag, that is, it may affect 
the boundary layer flow very slightly. 

The heat transfer on the nose of the airfoil is most 
easily comparable with that on the pipe in oblique flow, 
but it is just as difficult to treat theoretically as that 
on the front of a pipe. However, the following holds true: 

a) The heat dissipation of a surface area is so much 

greater as it is closer to the stagnation 
po int ; 

b) It varies in relation to the angle of attack and 

the Reynolds number similar to the dissipation 
of the side of the wing of which it occupies 
the major portion; 

c) The smaller the heat dissipating area near the 

stagnation point is, the greater its heat 
dissipation will be and so much less its de- 
pendence on the angle of attack. Pairly 
small areas fall within the laminar region 
"(approximately Rex<10 ); hence JTux varies 
verv slowly with the Reynolds number (with 
Rex 1/g ). 

With exception of the tests on the heat dissipation of 
wing radiators, described in R.&H. Report Fo. 1163, only 
mean values for the ITusselt number are known. But within 
the framework of the necessary experimental data the distri- 
bution of heat dissipation along the outer skin of the sur- 
face by constant limitation of the skin radiator is of par- 
ticular interest, because the local cooling conditions and 
the heat stresses are governed by it. If the analogy between 
flat plates and airfoils can be extended, the local values 
of heat-transfer coefficients can be closely approximated 
in power form STux = MRex^ from the measured mean values by 
differentiation of the empirical data, analogous to the func- 
tion represented in equation ( 19 ) , at least for the case 
where the heating extends without interruption from the lead- 
ing edge up to depth x. After that, it is important to 
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know- the appropriate areas on the wing surface for the 
heat dissipation in order to he ahle to find the "best 
arrangement for the particular case. 

For the heat dissipation of partial surfaces the 
statements made concerning the flat plate hold true. For 
every chosen parameter x x /x a factor w is computed 

similar to equation (22), with reference to the proper 
Nuxrp or Uux L ; w, which again can "be a functibn of 

Rex, might at the same time include the temperature ef- 
fect. Accordingly, it will he advisable to use the rela- 
tions governing the flat plate — extrapolated to increase 
of temperature 3 = 0 — as "standard" for a oasis in form 
of the curves Nux L , ETuxrp, and Nux{f. The actual vari- 
ation must, of course, he estahlished hy very careful 
experiments, while the concept "transitional region" re- 
quires a more exact interpretation. Comparison might also 
he made with a "standard function" for which ITuxjr is 

prohahly cost suitable , since it is free from empirically 
defined values. For reduction to other temperature ranges 
the previously introduced mean value for Pr should he 
retained and the variation of Pr in the expression V 
allowed f or . 

Irrespective of what new information is obtained for 
practical calculations the representation according to 
figure 1 should even then he more convenient than that 
according to figure 6, if the temperature effect is in- 
cluded in the suggested form. 

VII. SUPPLEMENTARY DRAG- OP SKIN RADIATORS 



The transition from the speed of undisturbed flow to 
zero value of the gas particles adhering to the surface 
takes place, as is known, within the boundary layer. In 
accord with measurements we chiefly differentiate between 
the parabolic velocity distribution in the laminar boundary 
layer and the distribution in the turbulent boundary layer 
according to the l/7 power law (pt. IV). If the tempera- 
ture of the surface differs from that of the gas, the 
statements concerning the temperature variation in the 
boundary layer hold true. If the relation of density p, 
kinematic viscosity v , and temperature conductivity a 
to air temperature is discounted, the fields of velocity- 
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and temperature are exactly similar for the specific case 
that friction and heat dissipation "begin simultaneously 
and the Prandtl number is equal to 1, the increase of 
temperature being computed from the surface temperature 
I 0 < Hany problems can be sat is factorily solved with the 
aid of this simple, although not exactly correct rela- 
tionship. For more accurate adaptation to practice the 
previously cited values \|/ and cp deduced from tests 
we used. 

The nature of the temperature effect on the drag con- 
sists, accordingly, in the simultaneous influence on the 
boundary layer thickness and further on the velocity pro- 
file within the boundary layer. In the test the two ef- 
fects can, of course, be measured only concurrently. 
Application to geometrically similar processes is very 
reliably possible if the cr i t i cal 'values are borne in mind. 
But it first needs to be proved that the data on the drag 
of thin, heated plates is applicable to airfoils; the ef- 
fect on the boundary layer thickness being treated as if 
it did not change the velocity distribution within the 
boundary layer. 

It may be suspected that with regard to drag a heated 
airfoil is decidedly thicker by £Ae than the cold one, 
and the conversion of the thereby increased flow veloci- 
ties to pressure could favor the formation of further form 
drag. Owing to the variation of the increases of velocity 
past the curved surface the boundary layer will form some- 
what different from that of a flat plate but the discrep- 
ancies are small. (Of. pt . V.) Undoubtedly the form drag 
(local vortices) is also of some influence. 

In order to make the estimation of the hypothetical 
thickening clear, an airfoil with vanishing form drag but 
finite skin friction is investigated, in its effect on a 
potential flow, first cold, then heated to temperature in- 
crease 0. 3y assumption of the decisive velocities v x for 
both surfaces of the wing the boundary layer thickness on 
all points of the airfoil circumference is (reference 2, 
pp. 138-139). 



5 



5.477 Hex' 



— 0.5 



-For laminar boundary layer 



(28) 



5 



0.37 Hex' 



— 0.2 



For turbulent boundary layer 



(29) 
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An. airfoil without skin friction — that is, on which 
the potential flow remained up t o the surface would, at 
equal pressure — and hence, Velocity distribution, have 
to he just as much thicker as the "boundary layer thickness 
decreased, if the potential velocity v x prevailed in it. 

For the thickness fl' of the substitute layer — if 



figure 8 is the distance of a filament from the 
within the boundary layer — must he 



y 

surface 



m 



p u 8' 



P v xy dy 



According to general law of velocity distribution 



xy 



m 



where m = l/2 for laminar, and m = l/7 for turbulent 
boundary layer , we get 



1 + m 



6 - 6 1 = = 



m 



1 + m 



> 



(30) 



The new profile obtained on this basis, which similar 
to the so— called half— body extends to infinity, can of 
course, never exactly fulfill the theoretical condition for 
zero skin friction (sum of all pressure forces equal zero) 
because of the already mentioned simplifications, and par- 
ticularly, because the potential velocity v x is not quite 
reached on the "outer limit" of the boundary layer. 



On the heated profile the boundary layer thickens up 
as a result of the heat expansion. Since no appreciable 
acceleration forces can occur in flow' direct ion by reason 
of the small inclination of the outside limitation of the 
boundary layer toward the surface, it may be assumed that 
each infinitely small width dy expands almost exactly 
in correspondence with its temperature. Hence, according 
to f i gur e 8 : 
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V \ 



m 



6 Q - 



« 5 



1 + m T 



m 0 



f 1 + — - ) 8 

V 1 + a I 



The substitute thickness 68 1 for the boundary layer con- 
verted for the potential velocity v x follows from 



P v 5 e' = J Py v xy Q-7 



T T 
P — v xy -=f dy 

xy X 



0 

= J* P v xy 



dy 



8q« = 



1 + in 



5 ( = 5') 



(31) 



To transform the actual profile in a "equivalent" 
profile without skin friction it must be thickened at all 
points by £9 or, compared to the cold one, by AG. 

According to equations (3l) and (30) and figure 8: 



m 



e 



AG = Go — g = — 

e* 1 + m T 



(32) 



For the general case of a profile with finite form drag 
this argument holds true at least in its tendency. Hence 
the following conclusion: for the" flow past the profile 
the latter on heating — as compared to cold — seems to be 
thicker by the sum of the quantities AG for the two 
profile surfaces. Hence, with equation (28): 
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EAS = S — — 5.477 Hex °" 5 x Tor laminar boundary layer (33) 
3 I 

with equation (29) 
ZAe = 2 -|- 0.37 Rex -0 ' 2 x For turbulent boundary layer (34) 



The inferior accuracy of equation (32) justifies, 
for simplification, the introduction of the undisturbed 
velocity v in place of the decisive velocities v even 
at greater angles of setting and a mean value instead of 
the arc lengths x computed for both wing surfaces from 
the stagnation point to the point of maximum profile 
thickness. Lastly, it is assumed that the distance of the 
point of maximum profile thickness is about 30 percent of 
the chord length t from the wing leading edge. Since 
the profile circumference is slightly greater than 2t 
it is approximately 

Hex = 0.15 Ee 
x = 0.3t 



The thickening XAe is referred to the actual maximum 
profile thickness d itself: 

Vd = ^AL 
d 

Then the equation for the laminar boundary layer, (32) reads: 
Y dL - 2 "J- ■?- 5.477 ( 0.15a©)" 



clL 



3 t 

~ 2.8 -L 
d 



°' 5 °' 3 4 

d 



Ee" 



■ 0.5 8 

I 



> (35) 



and for the turbulent boundary layer, equation (34): 
v dT ~ 2 T °-37 ( 0.15 Ee) -0,2 0.3 |- 



T d!r ~ 0.041-1- Ee" 0 - 3 f- 



(36) 
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Bquation (35) will serve for small models, equation (36) 
for full—scale versions. In all practical cases the 
amounts of 7& lie at 1 percent,* hence the form drag 
is scarcely affected by the heating in normal flight 
attitudes. Therefore the test data on flat plates can 
"be applied fairly reliably to airfoils. In this connec- 
tion it is further of interest to know whether portions 
of the heated outer skin lower the lift of the profile, 
especially its maximum, and whether the flow is adversely 
affected, so as to induce earlier "breakaway than on the 
unheated wing with increasing angle of setting. Such 
phenomena are said to have been observed in past experi- 
ments and they appear altogether plausible according to 
the foregoing arguments, even though these do no longer 
permit of safe conclusions in these extreme cases. If 
th* flow still adhers on the cold wing the assumed slight 
thickening up of the profile due to heating nay very well 
be sufficient to make it separate. 

Despite the fact that Ten Bosch's method appears to 
make the temperature effect excessive, the lack of other 
data prompts its use for determing the skin friction on a 
thin, flat plate heated over its total length to uniform 
increase of temperature P) and surface temperature T 0 
in two-dimensional flow. This merely involves the con- 
version of Hex to Rex 0 by means of equation (27) and in- 
sertion in Cf along with the still unknown factor \ , 
the- numerical values of which arc taken from pipe flow. 
The ratio of resistance of the heated to the" cold plate 



Cf-e 



/ „ \ — 0. 5 

for "short" plates (predominately laminar boundary layer, 
E.ex 0 <Sex cr ) according to equation (ll). 

b) T ctJ = 1 (38) 

for "medium" plates, (largely transitional boundary layer , 
He:c cr «Sex 0 <2X10 7 ) according to equation (14) 

*5 , or the small model cited in the subsequent footnote, we 
get in round numbers Re ~0.7 X10 6 (v~ 20 m/s, as minimum 
value), 8~75°, T~298°, d/t -0.13. According to equa- 
tion (35) V dL = ~0.01 (1 percent). If Re~10 7 , 6-100°, 
T~223°, d/t -0.13. Equation (33) gives 7 = -0.01 
( 1 per cent ) . ' 
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for "long" plates (predominately fully turbulent boundary 
layer, Rex » 2 X 10 7 ) 'according to equation (16). 



With the approximated relation 




these 



facts become even more evident. 

The evaluation of equations (37), (38), and (39) is 
graphically represented in figure 9. On the basis of the 
little reliable assumptions made, it is found that: on 
short plates the skin friction drag rises with ascending 
8 but slower for high temperatures, with an initial 
drag decrease at low 8; the same holds for long plates, 
but the initial decrease is more distinct, and the curve 
is subsequently flatter. On medium plates, on the other 
hand, the frictional drag drops to a minimum value, with - 
ascending 8, independent of the actual height of tempera— 
tur e . * 

If the fore part of a plate is not heated, its skin 
friction drtig naturally decreases proportionally less 
than corresponds to the ratio of the heat volumes (factor w). 
The exact magnitude of the frictional force exerted by the 
flow through the heated boundary layer onto a cold portion 
disposed behind the heated surface must be decided by ex- 
periment. Numerically this proportion of the drag lies 
naturally between that occurring in isothermic flow, if the 
temperature is, once T, then T Q . 

*The heat dissipation of surface portions distributed over 
the nose, upper and lower surface to the amount of ~55-*per— 
cent of the profile circumference was measured on a model 
wing, beatable to 250 millimeters chord (0. Seibert, Jahr— 
buch 1938 der deutschen Luf tf ahr t f or s chung , pp. II 224 — 
II 244). Subsequent drag measurements, by 0. Schoppe, 
disclosed as definite increase in the total drag on heat- 
ing throughout the entire possible range (Re to ~ 2 X 10 s ) . 

The average for Vc was 1.12 at to "100°; t ~ 25° 
(8-75°). The proportion of the laminar flow to the skin 
friction was certainly still perceptible on the model. 
Greater geometrical similarity would have perhaps afforded 
a lower Y c at these temperatures. (Continued on p. 30) 
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On a wing with, heating surfaces in the foreward portion 
only, the pressure distribution measurement in the rear 
unheated portion which can he easily made during the heat 
transfer experiments, itself is likely to afford practical 
data for the evaluation of the flow and separation phe- 
nomena. 

Impending measurements may yield smaller drag varia- 
tions than those deduced. But in any event the tenden- 
cies will hold good to some extent. Perliaps this dis— 
similarily in "behavior of plates or airfoils depending 
upon the momentary Reynolds number enables an explanation 
of several contradictory observations on the relationship 
between drag and temperature. 



Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 



*(Continued from p. 29) Figure 9 shows V c at about 1.06 

for the cited conditions. If the finite 6 were exactly 
allowed for, it would have had to result in a value higher 
than that recorded, since completely heated surface had 
been assumed for the study and the temperature effect was 



1.22 actually greater than the test value V 0 . This seems 

to suggest that the unknown factor i is considerably 
closer to 1 for flat plates and airfoils than for the flow 
in pipes. 



known to be too great. 'The term 
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Figure 1.- Representation of Nu » f (Re) according to R and u 
Reports 1163, 1311, 1326, and 1481. 
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Pigure 2.- | values for air. 



Figure 5.- 



v x /v (for cos 
3-1) total 
plate plotted 
against c a 
(according to 
equation 24) . 
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Evaluation of equation (19; 



Fig. 3 
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Subscripts: 
, L = laminar "boundary layer 
U s transitional region 
Pr = (according to Prandtl) 
T * turbulent boundary layer 
K = (according to Karman-Latzke) 

Figure 3. -Relation of drag coefficients Cf and the mean Nusselt 

number Nux for the flat plate with the Reynolds 
number (increase of temperature =. o). i 



NACA Technical Memorandum No. 1044 



Fig. 4 




Figure 4.- Comparison of the theoretical relation Nux and 

Rex of the flat plate with Ellas' test data 
xi depth of cold fore portion 

x depth of total plate 
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Fig. 6 
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Figure 6.— Reduced actual Nux plotted against the theoretical 
Kux for the flat plate. 
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Fig. 7 
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RAF 266 « 300 
RAF 306 = 73° 

Figure 7. Figure 6 rep lot ted against Pe 0 to allow for 

the temperature effect according to ten Bosch 
(surface temperature of flat plate the same as that of 
mode Is ) . 
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Figs. 8,9 
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Figure 8.- Effect of heating on the flow, 

6; 69 = .actual "boundary layer thickness 
6';5e'= substitute thickness 
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Figure 9.- Drag variation of plate due to heating = Vc 
' _ _ (assumptions) 



